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We propose a phenomenological model of vortex tube reconnection at high Reynolds numbers. The 
basic picture is that squeezed vortex lines, formed by stretching in the region of closest approach 
between filaments, interact like dipoles (monopole-antimonopole pairs) of a confining electrostatic 
theory. The probability of dipole creation is found from a canonical ensemble spanned by foldings 
of the vortex tubes, with temperature parameter estimated from the typical energy variation taking 
place in the reconnection process. Vortex line reshuffling by viscous diffusion is described in terms of 
directional transitions of the dipoles. The model is used to fit with reasonable accuracy experimental 
data established long ago on the symmetric collision of vortex rings. We also study along similar 
lines the asymmetric case, related to the reconnection of non-parallel vortex tubes. 



There is growing evidence that the dynamics of vortex 
tubes is a necessary ingredient for a deeper understanding 
of turbulence. This view, initially suggested by images 
of the vorticity field produced through direct numerical 
simulations [|lj, received strong support from the recent 
accurate determination of scaling exponents for the veloc- 
ity structure functions, within a phenomenological theory 
which places filamentary configurations on a fundamen- 
tal status iQ. However, the present knowledge on vortex 
dynamics is still far from being complete, so that even in 
simplified situations, as in the collision of vortex rings, 
a formal theory remains to be developed. This difficulty 
is related in part to the absence of comprehensive phe- 
nomenological descriptions that could provide a starting 
point for more elaborate discussions. 

We will focus our attention on the problem of vortex 
reconnection. Its importance - considered more as an 
expectation in the long run ~ relies on the idea that the 
global structure of turbulent flows may depend on topol- 
ogy changing processes, like the intertwining of closed 
vortex tubes. Previous theoretical studies on vortex 
reconnection |^-^ basically correspond to the case of 
low and moderate Reynolds numbers. These attempts 
may be regarded as the counterpart to interesting ex- 
perimental observations reported by different groups on 
the colhsion of vortex rings |@-^. On the other side, 
at higher Reynolds numbers, relevant effects come into 
play, as stretching and the possible existence of singular- 
ities . Another important element to be considered 
at high Reynolds numbers is that vortex tube evolution 
is hardly reproducible, due to sensitivity to initial con- 
ditions, and one has to resort, therefore, to statistical 
methods. While it could seem there is no hope of an 
analytical treatment, since there are no standard tech- 
niques to find statistical measures in unstable dynamical 
systems, we show in this note that the well-known Boltz- 
man distribution appears as a natural candidate from 
which plausible consequences may be derived. 

In an experiment performed about 25 years ago, Fohl 
and Turner [Q, studied the symmetric collision of two 
identical vortex rings in water at high Reynolds number 



{Re ^ 4000), approaching each other with variable an- 
gle 29. They found that the fusion of colliding vortex 
rings, both with radius r and velocities (0, v sin 9, v cos 9) 
and {0, —V sm9,v cos9), into a single ring with velocity 
^ (0,0, —w/2) and radius 2r allways occurs in a first 
stage. The fused ring exhibits amplitude oscillations 
|l2| , [l3| , so that a second stage characterized by a splitting 
reconnection takes place with probability p = p{9) . The 
two rings created after the second reconnection move in 
a plane perpendicular to the initial collision plane. An 
important feature in the experiment is the existence of a 
critical angle. For 9 > 9c — 16°, it holds p{9) = 1, where 
as p{9) — > as 6* — > 0. An explanation of 9c was given by 
Fohl and Turner, based on the fact that the modes of am- 
plitude oscillations, which describe perturbations around 
a vortex ring of radius r and velocity v, have wavelength 
A„ — 2'Kr/n and period 



Tn{r,v) 



2-Kr 



n{n^ — 1) 2 1; 
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with n > 2. In a collision defined by the angle 29, the 
projection of the ring's velocity on the direction trans- 
verse to the symmetry plane is vsm9. One may expect 
the amplitude of oscillations in the fused ring, which de- 
pends on the collision angle, to be 2r when 



2r 



V sm u = 
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In this case, diametrically opposite parts of the fused ring 
will touch each other, producing reconnection. Eq. (|^) 
may be readily solved, yielding 9c = arcsin(\/3/27r) ~ 
16°. 

There are important questions not answered on this 
problem. We would like to understand in a more detailed 
way the form of p{9) , taking into account its behavior at 
small angles. As we show below, this may be achieved in 
an effective model where an important role is played by 
the dynamics of interacting dipoles. 

To start, imagine two vortex tubes, locally antiparallel 
in some neighbourhood Q, both carrying the same flux 
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and having identical circular cross sections. A descrip- 
tion of the physical mechanism underlying reconnection 
was put forward by Saffman ||] . In his model, the strain 
field shrinks Q in the plane perpendicular to the vortex 
tubes, so that viscous annihilation of vortex lines occurs, 
reducing the circulation in fl. Therefore, a pressure gra- 
dient along the vortex's cores develops, which increases 
strain and then viscous diffusion, enhancing reconnection 
in a self-induced way. The equations assumed to describe 
these steps give reasonable answers, in spite of some dis- 
agreements with real and numerical experiments [^,^. 
Saffman's model is in fact devoted to the situation of 
two antiparallel vortex tubes interacting at close enough 
distance. The model works better in the case of strong 
viscous diffusion (low Reynolds numbers), where vortic- 
ity amplification is not very large. 

We suggest here a scenario of reconnection at high 
Reynolds numbers, when strecthing effects become rel- 
evant, which probably does not disagree with Saffman's 
model, since its application will be related to a differ- 
ent range of physical parameters. The picture we will 
consider is that in the collision of vortex tubes a sys- 
tem of "dipoles" , emerges after some stretching in a pro- 
cess characterized by very small energy lost. Reconnec- 
tion is finished with right-angle transitions and subse- 
quent collapse of dipoles, through vortex line reshuffling 
by viscous diffusion. The main events are depicted in 
Fig. 1. It is important to keep in mind that dipoles are 
just effective structures which represent stretched vortex 
tube segmentsl^l In a more rigorous approach the recon- 
nection problem should be addressed in terms of vortex 
sheets rather than quasi one-dimensional supports of vor- 
ticity, since numerical simulations P,p^ show that vortex 
tubes flatten in the reconnection region due to stretching. 
Dipoles have to be regarded only as an useful approxima- 
tion, from which it is possible to get phenomenological 
results in the simplest computational way. 

We are interested in studying the energy of a flow given 
by two long vortex tubes, both carrying vorticity flux 0, 
with stretched segments of length S and separated by 
the distance d, as shown in Fig. 1. Vorticity may be 
decomposed as = u)^^'' + uj^^\ where u)^^'' is the field 
locally amplified by stretching, and w*^^^ is the field at 
other places along the tubes. The energy E — E{d, d) 
may be written as 



E — El + E2 + E12 
where, introducing the notation 

1 /" .Q 1 
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we have 



*We mean, throughout this work, an analogy with the 
electric dipole definition. 



El = (u;(i),J«) , E2 = (cJ(2)^^(2)^ ^ 

^12 = 2(c3W,c3(2)) . (5) 

To simplify expressions, the fator I/Stt in (jj) will be 
supressed henceforth, corresponding to the replacement 
(f) 2\/27r</). Taking the stretched segments to be iden- 
tical, both with circular cross section of radius e <ti S, it 
follows that 

^=ln(-)-arcsmh(-)-l + [(-)^ + l]^--. 

(6) 

In the computation of Ei , the kernel in (^ is regularized 
by means of 
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\x-x'\ [{x-x'y + e^]i 



(7) 



The contribution to the energy which comes from the in- 
teraction between the stretched and non-stretched parts 
is 



£;i2=(j(i^),l3(2^)) + ((:;(i^),c3(2«)) 



(8) 



with the superscripts L and R denoting the left and right 
vortex tubes. The effect of screening between the two 
vortex tubes is evaluated from the following estimative 
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The above expression is derived from a discretized ver- 
sion of the integral (^) in terms of vortex tube segments 
of length (5, as discussed in ref. [|l5|, where we additionally 
used reflection symmetry (invariance of the vorticity fleld 
under the interchange L <-> R. In the reconnection exper- 
iments, S typically fluctuates around some mean value 
6 > d. This inequality implies, with (^ and (P), that for 
e/5 <C 1 the interaction energy E12 may be neglected, 
when compared to Ei, and, therefore, E ~ Ei + E2 (as 
an example, consider the choice e/5 ^ 10~^, which gives 
Ei2{S,d)/Ei{6,d) < 0.1). We conclude that the positive 
energy variation due to stretching is local, being com- 
pensated by reduction of energy through foldings of the 
vortex tubes |l^ . 

There is a close relationship between the stretched vor- 
tex tube segments and dipoles of a conflning electrostatic 
theory. Let pi (x) and P2 (x) be the dipole moments of two 
charge distributions, respectively defined in compact re- 
gions fii and Q2- The interaction energy associated to a 
linear confining potential is 



Eir, 



d^xd^x'\x - x'\\7 ■ pi{x)\7 ■ p2{x') 
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Through integration by parts we get 
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{pi,P2) - {Pi{x) ■ r)ip2{x') -rjlr]' 
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where r — x—x'. Taking fii and fl2 as oriented segments, 
parallel to the vorticity field and parametrized by (0, 0, s) 
and (— d, 0,6 — s), with < s < 6 and 6 <^ d, the scalar 
products in the second term in eq. ( pl| ) may be neglected. 
In this situation the result may be identified with the 
interaction energy associated to vorticity fields given by 
pi and p2 ■ 

An alternative and straightforward way to obtain the 
connection with confining electrostatics is to represent 
the dipoles as "monopole-antimonopole" pairs. This is 
done by replacing the stretched vortex tube segments by 
monopoles at positions (0, 0, S) and d, 0, 0) and anti- 
monopoles at (0,0,0) and {—d,0,5). These points are 
just the boundaries of fli and $12- The monopole (an- 
timonopole) is the source (sink) of a radially symmet- 
ric vorticity field, which is not solenoidal - and hence 
deprived of direct physical meaning. The field of a 
monopole-antimonopole pair is given as w = +uj 
where 



(-) 



X±\ 
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Above, x± gives the position of the monopole with charge 
±0. To compute the energy of a system of monopole- 
antimonopole pairs, it is necessary to regularize infrared 
divergencies. Defining the flow inside a large sphere of 
radius A — )■ oo, we will have 



) = 2?7(p,g)r(A- Ifp-x^l) 



(13) 



where p,q — ±, and 77(73, q) = I for p = q, otherwise 
r]{p,q) = —1. The energy of the interacting monopole- 
antimonopole pairs considered here is, thenH, 



(14) 



When 6 ^ d, the above expression differs from eq. 
(^, only by self-energy terms, which are independent 
of d (for d/6 > 1 the agreement is better than 90%). 
Since the interaction potential is linear, the force between 
monopoles has constant strength and is directed along 
the line joining the particles. This suggests, in the de- 
scription of stretched vortex tube segments as monopole- 
antimonopole pairs, that whenever 6 > d, reconnection 
takes place. This is in fact the same condition that fol- 
lows from El (S, d) > Ei (d, d) , using the more precise re- 
sult, eq. (^. These are the energies for the configurations 
shown in Fig. 1. Furthermore, the energy dissipated in 
the reconnection process is AE = Ei (5, d) — Ei {d, 6) . 

To model the symmetric collision of vortex rings, it is 
necessary to know how d depends on the collision angle. 



^Note that the infrared divergencies cancel for a neutral sys- 
tem of monopoles. 



This may be obtained by replacing the numerator 2r in 
eq. (H) by 2r — d/(2ci), which we propose to be the oscil- 
lation amplitude at angle 9, where ci is a phenomenolog- 
ical parameter. We find d{9) = 4cir(l — 27rsin6'/V3). As 
only small angles are involved {0 < 9 < Oc — 16°), this 
expression may be effectively thought as a linear interpo- 
lation between the maximum value d = 4cir, which oc- 
curs at 9 = 0° , and the minimum value d = 0, at 9 = 9c- 
Another important parameter is the mean length of the 
stretched vortex tube segments. We define it as 5 = C2r. 
It is convenient to use dimensionless units where 6—1 
and 



d(0) ^4^(1-^ sin ( 
C2 V3 



(15) 



Due to the additive property of energies, E = Ei + E2, 
we may interpret the random behavior of 6 as being re- 
lated to fluctuations of Ei in a canonical ensemble. The 
elements of the canonical ensemble correspond to folded 
configurations of the non-stretched parts of the vortex 
tubes, which act like a reservoir exchanging energy with 
the stretched region. The probability density to have 
stretching length < 5 < 00 is 



p{S,d) ^ Z-^exp[-PEi{S,d)] , 
where f3 is the "inverse temperature" and 

Z = I dxexp[-(3Ei{x,d)] . 
la 



(16) 
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is the partition function. The definition 5 = 1 is used 
to find the temperature, /3^^ ~ Ei{1.0,d), which is the 
energy necessary for the creation of stretched segments of 
length 6. The probability to have reconnection is, thus. 



p{9) = / dxp{x,d{e)) 



d{9) 



(18) 



It is not our purpose to derive the above statistical me- 
chanical correspondence from first principles; we take it 
as a hypothesis. The main problem would be to show 
that energy fluctuations take place in a time scale much 
larger than the one of reconnection /4>, see Q). On 
the other hand, the relaxation time for the vortex sys- 
tem to reach thermal equilibrium has to be much smaller 
than the time spent for the whole process (~ r'^ /(f), last- 
ing from the fusion of the vortex rings up to the instant 
of split. Amplitude oscillations may be an important as- 
pect of such an analysis: at large n, the wave velocity is 
An/Tn ~ nv, according to eq. (|^). Therefore, "thermal 
equilibrium" would be assured by the fast propagation 
of perturbations along the vortex tubes. Actually, this 
equilibrium is not stable, due to the attraction between 
dipoles. It is worth noting that a more complete study of 
the right-angle transition of dipoles would probably deal 
with methods of non-equilibrium statistical mechanics. 



3 



We need to set the value of the dimensionless phe- 
nomenological parameter C\jci. This is done in princi- 
ple by searching for the best agreement with experimen- 
tal data, but we do not expect this parameter to sig- 
nificantly depart from unit. The reason is that 5 ~ r, 
as indicated in numerical and real experiments, and 
(i(0°) ~ 4r, if one assumes that the amplitude of os- 
cillations in the fused ring vanishes when 6* = 0°. We 
note that d(0°)/(5 — ^cijci is a quantity suitable of ex- 
perimental determination. The resulting is shown 
in Fig. 2, with Cxjc^ = 1.0 and e = 0.01. As a matter 
of fact, the form of p{d) is not altered in an important 
way for e < 0.1. In the limit e/6 0, we may solve 
eq. ( p^ exactly, to get that the probability of reconnec- 
tion decreases exponentially with the distance between 
the vortex tubes, that is, p{6) = exp[— 

We may proceed through similar computations to 
study the asymmetric scattering of vortex rings, aiming 
at predictions that could be tested in future experiments. 

In the symmetric case considered so far, let A be the 
point where the colliding vortex rings first touch. A is di- 
ametrically opposite to points B and C, which are placed 
in different rings. The angle at the vertex A, defined by 
the segments AB and AC is 180° - 29. We study the 
asymmetric collision obtained from the following initial 
configuration: while the ring which contains C is fixed, 
the other ring is rotated around the axis AB by the angle 
a. This is precisely the angle between the vortex tubes 
at the point of contact. 

We want to find now the probability p{9, a) for the 
splitting reconnection. As the fused ring evolves, points 
B and C move toward each other with relative velocity 



(19) 



2usin6'cos^(a/2) 



This amounts to replacing sin 9 which appears in eq. (|1E^ 
by sin6'cos^(a/2), to define d{9, a), the distance between 
the non-parallel dipoles. Assuming that the second re- 
connection occurs close to points B and C, with stretched 
vortex tubes keeping the initial relative angle a, all we 
need to know is the energy of such a configuration. Ne- 
glecting the interaction terms in the expression for the en- 
ergy, which depend on the distance between dipoles, the 
condition for reconnection is then Ei{5) > Ei{t{9,a)), 
where 



9, a) = [£(9,a)+8'^ sin2(a/2)]- 



(20) 



Therefore, reconnection occurs when 5 > l{9,a), or, 
equivalently, 5 > d{9,a)/ coa{a/2). We just have to 
replace the lower bound d{9) in the integral ([l^ ) by 
d{9, a)/ cos(a/2). We carried out computations using the 
same set of parameters (3 and ci/c2 for the former case 
{a — 0). It is possible that deviations grow as a gets 
larger, where the dipole model may loose its applicabil- 
ity. However, there is some indication that reconnection 
is supressed at large a , which is also verified through 
an explicit computation of p{9,a). In Fig. 3 p{9,a) is 
shown with < a < 180°, for 9 = 15° and 9 = 17°. The 



latter situation is perhaps the more interesting, because 
of the plateau given by p{9^ a) = 1 at small values of a. 

To summarize, we studied the problem of vortex re- 
connection at high Reynolds numbers, where stretching 
effects become important. A simple correspondence with 
confining eletrostatics and statistical mechanics allowed 
us to investigate the way how "dipoles", i.e., stretched 
vortex tube segments, behave in the process of reconnec- 
tion. The initial dipole configuration evolves, thanks to 
diffusion, towards a state which links the interacting vor- 
tex tubes. The measurements of Fohl and Turner Q for 
the probability of a splitting reconnection after the ini- 
tial merger of vortex rings are successfully accounted by 
the present model. We also defined some predictions con- 
cerning the case of asymmetric collisions, to be compared 
with possible future experimental observations. 
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FIGURE CAPTIONS 



Fig. 1 

Sketch of the configuration of the vortex tubes in 
the region of closest approximation. An attractive 
force betwen the oppositely oriented stretched segments 
( "dipoles" ) is followed by vorticity cancelation (not rep- 
resented in the picture) and then by the right-angle tran- 
sition of the dipoles. 

Fig. 2 

The probability p{9) for the occurrence of a splitting 
reconnection. Dots represent the values measured by 
Fohl and Turner. The continuous lino is the predic- 
tion of the dipole model obtained with /3 = 1/ Ei{l.Q,d), 
ci/c2 = 1.0, and e = 0.01. 

Fig. 3 

The probability p{9, a) for the ocurrence of a splitting 
reconnection in the asymmetric case, where 6 is fixed and 
< a < 180°. a) 6* = 15°. b) 6* = 17°. 
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